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Abstract 



The density functional determining the Coulomb energy of nuclei is calculated to 
the first order in e^. It is shown that the Coulomb energy includes three terms: the 
Hartree energy; the Fock energy; and the correlation Coulomb energy (CCE), which 
contributes considerably to the surface energy, the mass difference between mirror 
nuclei, and the single-particle spectrum. A CCE-based mechanism of a systematic 
shift of the single-particle spectrum is proposed. A dominant contribution to the 
CCE is shown to come from the surface region of nuclei. The CCE effect on the 
calculated proton drip line is examined, and the maximum charge Z of nuclei near 
this line is found to decrease by 2 or 3 units. The effect of Coulomb interaction on 
the effective proton mass is analyzed. 
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I. INTRODUCTION 



Relatively weak Coulomb interaction substantially affects the properties of nuclei owing to its 
long-range character. It is reliably established that the electrostatic Coulomb energy of the dis- 
tributed charge is the dominant term in the Coulomb energy of nuclei. This term, also known as the 
Hartree energy, is proportional to Z^e^/i?, where Z is the number of intranuclear protons, R is the 
nuclear radius, and e is the proton charge; at large values of Z, it induces the breakup of nuclei. At 
the same time, there are a few less significant contributions to the Coulomb energy. These include the 
Fock exchange contribution and a number of terms associated with the correlated motion of nucleons. 
It will be shown below that, in constructing the density functional determining the Coulomb energy 
of a nucleus, it is sufficient to calculate all these terms to the first order in e^. The presence of these 
terms is clearly illustrated by the well-known Nolen-Schiffer anomaly [|I|], which was discovered more 
than 30 years ago. This anomaly was deduced from a comparison of the calculated mass difference 
between two mirror nuclei with its experimental value. The comparison revealed that, on average, 
the results of the calculations performed without taking consistently into account the correlation 
Coulomb energy (CCE) fall short of relevant experimental values by 10% This discrepancy can 
be removed by introducing charge-dependent forces whose strength is assumed to be preset by 
this discrepancy. However, this is possible if we are sure that nuclear dynamics is taken properly 
into account and that all the Coulomb contributions are consistently included. A new mechanism 
that enhances the contribution of Coulomb interaction to the energy of the nuclear ground state and 
which is caused by the presence of a surface (a general property of equilibrium finite Fermi systems) 
was found in [§,0. This contribution, which is proportional to the nuclear surface (Z^/'^), made it 
possible to explain a dominant part of the above anomaly. 

The main objective of this study is to calculate consistently the density functional determining 
the Coulomb energy. As the result of this calculation, which will be performed to the first order in 
and which will take into account the mechanism enhancing the CCE contribution, we will be able to 
clarify the effect of this mechanism on the single-particle spectrum of nuclei, on the effective nucleon 
mass, and on the position of the proton drip line. In what follows, we do not highlight the difference 
between this functional and the Coulomb energy unless this leads to confusion. 

The ensuing exposition is organized as follows. A general formulation of our approach to calcu- 
lating the Coulomb energy is given in Section 2. Section 3 is devoted to calculating the systematic 
CCE-induced shifts of the single-particle spectrum and the proton drip line. The effect of Coulomb 
interaction on the effective proton mass is considered in Section 4. Section 5 summarizes the basic 
results. 

II. COULOMB CORRELATION ENERGY 

In order to formulate the aforementioned mechanism and to consider its effect, we will make 
use of the Hartree-Fock method employing effective forces, which was substantiated within density- 
functional theory 0]. In this case, the energy E of the nuclear ground state has the form 

E = Fo[pp, pn] + F^[pp, p„]. (1) 

where Fq is that part of the functional which is due to nuclear forces. Since we assume that the charge 
symmetry of nuclear forces takes place, Fq depends identically on the single-particle densities pp and 
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Pn of protons and neutrons (these densities are determined by minimizing E). The term Fc is due to 
other interactions, which are weaker than nuclear forces. The smallness of these interactions can be 
characterized by the ratio a of the potential energy of these interactions between two nucleons to the 
Fermi energy ~ 40 MeV. For the Coulomb interaction, this ratio is a = jr^Ep ~ 0.03, where 
ro is the mean spacing between nucleons in equilibrium nuclear matter. Therefore, it is sufficient to 
calculate the functional Fc in the first order in the relevant coupling constant. In the case of Coulomb 
interaction, which is the main subject of the present study, the functional F^ will be calculated to 
the first order in 62. This functional could be approximated by some simple expression, as was done, 
for example, for the functional Fq. However, experience gained in calculating the mass differences 
between mirror nuclei (these mass differences are directly determined by the functional F^ proves the 
inefficiency of such attempts. Moreover, this oversimplified approximation does not reveal the physics 
behind the functional F^.. Therefore, it is preferable to calculate the functional Ff.. In doing this, 
it is assumed that the functional Fq is known. It can be taken in the form specified by the Skyrme 
interaction or in the form proposed in |^ . The latter is advantageous in that it is characterized by 
a high precision, is simple in use, and involves the bare nucleon mass. By construction, the functional 
Fq must contain the unnrenormalized mass M of the proton (neutron) rather than its effective mass 
M* because the effective mass, as well as the nuclear spectrum, cannot be directly included in it 
rO|] . The single-particle distributions of the proton and neutron densities (pp and p„, respectively) 



are specified as 



Pp(r) 
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Pn(r) 
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where rip and n„ 
tively, while 0' 



' are the occupation numbers for single-particle proton and neutron levels, respec- 
and 0^ are, respectively, the proton and neutron single-particle wave functions de- 



termined from the Hartree-Fock equation §]. The term -Fc[Pp] is usually taken in the first order in 
the Coulomb interaction in the Hartree-Fock approximation 



^ 2 
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(3) 



Here, Xp(i'i) ^'25 1^) is the linear response function for noninteracting protons moving in the self- 
consistent single-particle field Vp. The first and the second term on the right-hand side of Eq. (3) 
are, respectively, the Hartree and the Fock density functional {F^ and F^). The exchange term is 
usually approximated by the expression 



3e2 /3\i/3 
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(4) 



However, the equality in (3) is not exact even in the first order in the Coulomb interaction (e^ order) 
because it does not include the functional F™'"^[pp(r)] taking into account the contribution to the 
Coulomb energy from the correlated motion of protons under the effect of the effective (residual) 
nuclear interaction i?a/3(ri, i'2)- In the first order in e^, the CCE functional F^"^"^ has the form 



Xpp(ri, r2, luj) - x^(ri, ra, luj) 



dridv2dui 
2^ ' 



(5) 
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where Xpp(ri,r2,w) is the exact function representing the hnear response of intranuclear protons to 
an external electric field. It should be noted that the function Xpp is completely determined by the 
functional Fq and that it involves no smallness associated with the Coulomb interaction. The 
total response function contains three independent components — Xpp? Xnn and Xnp — and satisfies 
the equation 



Xlm — X^l^lm + XI X^RlkXkr. 

k 



(6) 



where integration signs are omitted for the sake of simplicity. The response function fo^' nonin- 
teracting nucleons (protons) moving in the self-consistent field Vp has the form 



lr2j 
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The single-particle eigenvalues the frequencies u, 



ik 



e'^ — £^ 

^p ^p1 



the single-particle wave func- 



tions (/>p*(r), and the corresponding densities Pp and p„ are determined by solving the single-particle 



Hartree-Fock-like equations, or the Kohn-Sham equations fTT 




r), 



(8) 



which are derived by means of a conventional procedure that amounts to varying Eq. (1) with respect 
to the single-particle functions 0^ ^ [^|[ . The self-consistent potential Vp and the effective interaction 
Rai3 are given by 

/ N 5{Fo[pp,Pn] -To[pp,Pn]) 6^{Fo[pp,Pn] - To[pp, Pn]) , . 

= r-r-\ ! Raf3[ri,r2) = — T— ^ — , (9j 



6pp{r) 



6pa{ri)6pfs{ri] 



where To[pp, pn] is the kinetic-energy functional for noninteracting nucleons. Thus, we conclude that, 
to the first order in e^, the eventual form of the functional Fc is 



F, = + F; + Ff 



(10) 



Equations (1), (3), (5), and (10) provide a basis for consistently calculating the Coulomb energy and 
for taking into account the Coulomb interaction effect on the properties of nuclei. The poles of the 
response function c Xpp determine the collective spectrum of nuclei, while the residues at these poles 
govern the relevant transition probabilities. As can be seen from Eq. (5), a dominant contribution 
to F^°^^ comes from collective isoscalar surface vibrations whose excitation energies are much lower 
than those of corresponding isovector modes. Thus, the functional F^:"'^^ is controlled by the isoscalar 
components of the effective interaction. 

It is instructive to consider the CCE for symmetric infinite nuclear matter |^,^, in which case 
the relevant equations are considerably simplified. The function c Xpp assumes the form 



XppiQ,^^) 
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and determines the functional 



d(\duj 



(12) 



For symmetric matter, we have Xp = 
interaction /(p) between nucleons that was introduced by Migdal |]1T 



X°/2; the effective interaction R{p) is similar to the local 
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Considering that the isoscalar amplitude /(p 0) is approximately equal to —2.5 and that 
PfM/tt'^ = x°(0,zO), one can see from Eq. (11) that, at relatively low densities corresponding 
to the surface region of nuclei, the denominator 1 — -R+x° vanishes at small values of q and u. There- 
fore, the function Xpp develops poles on the imaginary axis. These poles, which evince instability 
of low-density nuclear matter |[Tl|, lead to the divergence of the integral in Eq. (12), and this pre- 
vents a calculation of the CCE at these nuclear-matter densities. We can conclude that the main 
contribution to this integral comes from the isoscalar response function. From the above, it also 
follows that the conventional procedure based on calculating various quantities in the local density 
approximation (LDA) in infinite uniform matter fails in this case. In this connection, it is worth 



noting that, although an attempt was made in |T2[ to construct the functional Fq within the LDA, 
this is impossible for the same reason. It is convenient to calculate the CCE in semi-infinite nuclear 
matter, where there is a region of low-density nuclear matter near the surface. In the generalized 
LDA, we then obtain 



[Pp 



Pp(r)ec(p(r))c/r, 



(14) 



where ec(p) is the CCE per proton and pp is the single-particle proton density. The calculation 
of the energy ec(pp) for semi-infinite nuclear matter revealed that, at the surface, this energy has a 
pronounced positive peak, which corresponds to smoothing the divergence in (12) for uniform nuclear 
matter M. It is convenient to approximate the energy Cc by the simple expression 



ec(pp(r)) = D 



Po 



Po 



4/3 



(15) 



where 2po = 0.16 /m~'^ is the equilibrium nuclear density and D ~ 6 MeV. In contrast to the 
terms F^ and Ff, the energy edpp), which has a pronounced peak at the nuclear surface, makes a 
noticeable contribution to the surface tension cr^, 



pp{z)[ecipp{z)) - ec{pp{-oo))]dz. 



(16) 



Hence, the Weizsocker mass formula must be supplemented with the term AE whose contribution 
to the total binding energy of a nucleus can be represented as 



AE = a[Z + a'^Z^/^, 



(17) 
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where a\ ^ —0.1 MeV and a'2 — 1.0 MeV. Thus, we can see from Eq. (17) that the surface tension 
has an isovector component. It will be shown below that this is so for the effective mass as well. 
It follows from Eq. (17) that the surface tension of the proton nuclear Fermi liquid is effectively 
greater than that of neutron matter. This means that, when protons are added to a nucleus, their 
binding energy becomes somewhat smaller because of an increase in the surface of the proton liquid 
component. It will also be proven below that this decrease in the proton binding energy is sufficient 
for explaining the mass difference between mirror nuclei. Further, the additional surface tension 
(Tc and the addition of extra neutrons are expected to induce opposite modifications in the root- 
mean-square radius < Vp > oi the intranuclear-proton distribution: the former would moderate an 
increase associated with the latter. The preliminary calculations reveal that this could explain the 
anomalously small increase in the radius < > in going over from the '*°Ca to the ^^Ca nucleus. 
The contribution (Xc must also be taken into account in considering fission barriers for heavy nuclei. 



III. SINGLE-PARTICLE SPECTRUM AND THE PROTON DRIP LINE 



Let us calculate the CCE-induced shift of the single-particle proton excitation spectrum . For 



this, we will use the well-known Landau equation 

5E _ I 



It follows from Eqs. (5) and (18) that this shift Ae' can be represented as 
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(19) 



The variational derivative 5x^/5^^ has the simple functional form 
'^Xp(ri,r2,cj) 



5nL 



GP(ri, r2, w + 4) + GP(ri, rs, -uj + e 



^(ri)0Ur2) 



(20) 



where G^{yi, r2, uj) is the single-particle Green's function for the system of Z noninteracting protons 
moving in a self-consistent single-particle nuclear potential. The functional derivative 5xpp/5n^p is 
determined by the matrix equation 



6nl 6nl 



Sim + 
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onl oni 



Xkm + Xl 
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(21) 



which can be derived by directly varying Eq. (6). Integration with respect to spatial coordinates is 
implied in Eq. (21) in just the same way as in Eq. (6). As can be seen from Eq. (9), the effective 
interaction i?^^ is determined by the form of the functional Fq. In order to simplify the calculations, 
we took, however, the interaction Rai3 in the separable representation [n 



Rikiri,r2) = X- 



dr 



dr 



(22) 
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where Vi{r) is the self-consistent single-particle proton (neutron) potential. The value of the param- 
eter A is chosen in such a way that the dipole response has a pole at = 0. The above separable 
form of the effective (residual) interaction is extensively used and, as was shown in [l7|,[T5|, provides 
a good description of collective nuclear excitations. The value calculated here for the shift ASp of 
single-particle proton levels located near the Fermi level is (0.2 — 0.4) MeV both for medium-mass 
and heavy nuclei, in overall agreement with the value of the Nolen-Schiffer anomaly. It is worthwhile 
to verify these results by using the simple LDA expressions (14) and (15) for The shift ASp 

then assumes the form 

^4= /^^?^l</>pWrrfr. (23) 



6p, 

It is convenient to approximate the single-particle density pp{r) by the Fermi distribution 

Pp^""^ ^) = TT m/ V ^24) 

1 + exp((r — it)/a) 

where R is the nuclear radius and the diffuseness parameter a is taken to be 0.6 fm. In this case, 
the functional F^^°^''"[pp] given by Eq. (14) reduces to a function F^°^^{R) of the nuclear radius. 
An addition of a proton increases the radius R by the quantity AR that is determined from the 
normalization condition 

I [ppir, R + AR)- pp{r, R)] dr = 1. 

The CCE-induced shift of a proton level occurring near the Fermi surface is given by 

Aep = F^"''{R + AR) - F™""(i?). 

After simple transformations of the relevant integrals of the Fermi functions, we find that, for medium- 
mass and heavy nuclei, the sought shift of a proton level near the Fermi level is 

Aep ~ ^ ^ (0.3 - 0.4)MeV. (25) 
2R 

Since corrections of order D{a/R)^ were discarded, the result in (25) slightly overestimates the shift. 

The calculated mass differences between mirror nuclei are quoted in the table. In this calculation, 
the functional Fq was constructed on the basis of the 5"/// interaction |^, while the functionals Fc 
and F^""^^ were, respectively, taken in the form (10) and specified by Eqs. (14) and (15). 



Mass differences between mirror nuclei 



(in MeV) 



Nuclei 


150 - i^N 


17F _ 170 


39Ca - 


41 Sc - 41 Ca 


48Ni - 48Ca 


Theory 


3.48 


3.56 


7.23 


7.24 


66.70 


Experiment |^,16|| 


3.54 


3.54 


7.30 


7.28 


67.06 



The table demonstrates that the mass differences between mirror nuclei are closely reproduced. 
The remaining disagreement can be used to determine the coupling constants for forces violating 
the charge symmetry of nuclear interaction. Thus, traditional nuclear physics can be advantageous 
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in calculating the constants concerning elementary-particle physics. Equation (25) shows that the 
shift Asp of a single-particle level approximately compensates for the level shift caused by exchange 
Coulomb interaction. It should be recalled that the exchange Coulomb interaction reduces the energy 
of a single-particle level, as can be seen from Eq. (4). Therefore, the total shift of a single-particle 
level due to Coulomb interaction can be derived by taking into account solely the direct Coulomb 
interaction — that is, by retaining only the term in Eq. (10). A similar procedure was postulated 
in iPJl^. It should be noted that the recipe assuming mutual cancellation of the exchange Coulomb 
interaction and the CCE does not lead to a noticeable value of the coefficient a'g in Eq. (17) (recall 
this coefficient determines the isovector component of the surface tension) because the contribution 
from and F^ [see Eq. (4)] to the surface tension is small. On the contrary, our calculations 
indicate that the coefficient 03 is not small. It also follows from Eq. (25) that, upon taking the 
CCE into account, the last filled single-particle proton level is shifted upward by 0.3 to 0.4 MeV. As 
a result, the maximal charge Z of a nucleus occurring near the proton drip line decreases by 2 or 
3 units, as follows from |17]. Therefore, we conclude that, in calculating the proton drip line, it is 
important to take the CCE into account. 



IV. EFFECTIVE MASS 



Let us consider the change AM that Coulomb interaction induces in the effective proton mass 
M*. In the case of uniform nuclear matter, the single-particle spectrum depends on the momentum 
p, while the effective mass, as follows from Eq. (18), is determined by the expression 



1 



1 dep{p ) 
Pf dp 



\p=PF1 



(26) 



where pp is the Fermi momentum. In order to calculate AM, we use Eq. (26), replacing ep{p) by 
the shift of the single-particle spectrum due to Coulomb interaction. We have 



AM 



d 6 



M* {M* + AM) Pf dp dup 



Xpp{q,iuj) + 27ippS{uj) 



dqduo 



|3' 



(27) 



Evaluating variations in Eq. (27) and considering that, in the case of uniform matter, Eqs. (6) and 
(21) reduce to algebraic equations, we obtain 0,|l^, 



AM 



e^d fSxoiq,iuj) 



dqdu 



M*{M* + AM) ppdpJ 6np {I - R+{q,iuj)xo{q,i^)f q^i'^T^f 



(28) 



where M* is the effective proton mass in the absence of Coulomb interact ion and Rj^ — (^Rpp-\-Rpn) /2 
is the effective interaction. The derivative d/dp was taken at the point p = pf- Let us consider the 
effective-mass variation AM for the case where the system under consideration is characterized by 
the compressibility parameter k, 1/k = p^d'^E / dp^ , tending to infinity. This state is similar to 
states existing in the surface region of nuclei. This consideration will enable us to draw qualitative 
conclusions on the behavior of the effective mass in a nucleus. It can be shown that UTI 



d_J_ 

dp 6n, 



'Xo{q,^)\ 



P^Pf 



T^iPF 

Pf 



|p + q|)5(^)p(p + q). 



(29) 
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Substituting (29) into (28) and retaining only the leading term, we obtain 



1 xdx , , 



M* + AM M* 27TPF J-i {1 -x)[l - R{q{x),0)xo{q{x),^ 



where q{x) = pi?y2(l — x). At the point where the compressibility diverges, the denominator (1 — 
R+XoY of the integrand in Eq. (30) vanishes at x = 1 (g = 0). The relevant integral, which is 
positive, diverges, which results in the vanishing of the effective mass — that is, M* + AM —>■ 0. This 
result demonstrates that, in equilibrium nuclear systems, the CCE noticeably affects the effective 
mass, the decrease AM in the effective proton mass being dependent on the degree to which the 
density of the system,p, deviates from its equilibrium value po- In a nucleus, the divergence of the 
integral in Eq. (30) — it could be due to the decrease in the density near the surface of the nucleus 
— gives way to a finite value. As a result, the effective proton mass decreases in relation to M* 
and becomes smaller than the effective neutron mass. It can be said that the effective nucleon mass 
acquires an isovector component |T^. This possibility was discussed in A theoretical validation 
of this effect has been given here. 



V. CONCLUSION 



The basic results of the present study can be summarized as follows. A consistent scheme for 
constructing the density functional determining the Coulomb energy of a nucleus has been developed 
to the first order in e^. Using this functional, we have calculated the single-particle spectra of nuclei 
and the systematic shift of these spectra that is induced by the CCE. It has been shown that the 
Nolen-Schiffer anomaly is removed to a considerable extent by this systematic shift of the single- 
particle spectrum. Owing to the same mechanism, the proton drip line undergoes a shift of 2 or 3 units 
toward smaller values of the charge Z of a nucleus occurring near this line. It has been shown that the 
CCE must be taken into account in calculating the effective nucleon mass. This opens the possibility 
of estimating the coupling constants for forces violating charge symmetry. The contribution of these 
forces must be treated in the same manner as this has been done for the Coulomb interaction. It 
can be expected that the proposed procedure will make it possible to construct a density functional 
applicable to describing various properties of nuclei occurring both in the valley of stability and 
beyond it. 
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